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This paper presents a new scheme for variable structure (VS) fuzzy PD 
controller. The rule base of the fuzzy PD controller is tuned online. The 
purpose of the proposed controller is to track accurately a preselected 
position command for the servomechanism system. Therefore, this study 
establishes a model using a black-box modeling approach; simulations were 
performed based on real-time data collected by LabVIEW and processed 
using MATLAB. The input signal for the servomechanism driver is a 
pseudo-random binary sequence that considers violent excitation in the 
frequency interval. The candidate models were obtained using linear least 
squares, nonlinear least squares, and deep neural network (DNN). The 
validation results proved that the identified model based on DNN has the 
smallest mean square errors. Then, the DNN identified model was used to 
design the proposed control techniques. A comparison had been executed 
between the VS fuzzy PD control, the conventional PD control, and the 


Variable structure fixed structure fuzzy PD control. The experimental results confirm the 
proposed VS fuzzy PD control can absorb the nonlinear behavior of the 
system. The speed regulation test, it reduces the rise time from 50% to 56%. 
While continuously changing in speed, it has the smallest tracking error 


(0.412 inches). 
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1. INTRODUCTION 

The one-stage servomechanism system is the basic unit of most CNC machines [1]. These types of 
machines need accurate motion control to produce high-quality products [2], [3]. Also, the high speed and 
high accuracy position applications have been rapidly and broadly developed to facilitate the complex 
automatic processes in the industrial field [4], [5]. Generally, the servomechanism systems don’t have an 
accurate estimated parameter for friction and backlash models which leads to system uncertainty. The friction 
causes high steady-state error while the backlash reasons undesired vibration in the system [4]. The models 
that can describe the friction phenomenon in servomechanism systems suffer from parameter uncertainty [6]. 
So, the control task for this operation faces a big challenge [7]. There are many design schemes for single- 
axis controllers such as conventional PID control, self-tuning proportional—integral—derivative (PID) control, 
fuzzy PID control, adaptive control, and sliding mode control [8]. The advantage of self-tuning and sliding 
mode techniques is that the tracking performance is good under the different operating points of position 
commands [9], [10]. In several applications, the sliding mode control has been applied such as microgripper 
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position/force control as presented in [11] and the position control of X-Y stage servomechanism as 
demonstrated in [12]. 

On the other side, the disadvantages are that the chattering phenomenon can be an obstacle to 
improve the control accuracy, and in some cases, it causes instability in the system [13], [14]. The fuzzy logic 
system is a suitable technique to treat the uncertainty and nonlinear problems of complex systems [15], [16]. 
Sometimes, the fixed rule base structure of the fuzzy system cannot deal with system uncertainty and 
nonlinearity so, several previous studies had attempted to use different techniques to develop an adaptive 
mechanism that can tune online the rule base or normalizing gains as in [17], [18]. However, the 
disadvantage of these mechanisms takes a long time to adjust accurately the rule base and need high 
processing, so it cannot execute practically [19], [20]. 

It is known that the model reference adaptive control (MRAC) is highly efficient adaptive control 
mechanism where it forces the overall system to follow the behavior of a designed model reference [21]. The 
selected model reference can be a first or second-order system according to the point of view of the designer 
[22], [23]. Experimental identification is a well-recognized methodology to obtain a precise process model, 
often intended for control but also other purposes. The nature of the input signal for the identification has a 
great effect on the accuracy of the model. But, in many cases, the input signal cannot be easily selected 
concerning plant behavior constrictions. pseudo-random binary sequences (PRBS) are often used as violent 
excitation signals for system identification, due to it has a finite length and can be synthesized frequently 
with simple generators while presenting favorable spectra for identification [24]. 

The family of candidate models for system identification can be classified along with several 
different aspects. The linear and nonlinear theory has been well developed and investigated for real 
applications throughout recent years. The linear structure is used to simplify the analysis where the 
parameters are constant and do not vary throughout a simulation, such as the autoregressive with exogenous 
inputs (ARX) model. In contrast, a non-linear model presents dependent parameters that are permitted to vary 
throughout a simulation run, and its use becomes necessary where interdependencies between parameters 
cannot be considered insignificant [25]. 

Lately, deep learning has become attractive and has significant attention from a wide range of 
engineering applications. Compare to traditional neural networks, the vital features of deep learning are to 
have more hidden layers and neurons and to improve learning performance [26]. Using these features, 
complex and large problems that could not be solved with traditional neural networks can be resolved by 
deep learning algorithms [27]. Therefore, deep learning has been subjected to various applications including 
pattern recognition and classification problems; for example, handwritten digit recognition, speech 
recognition, human action recognition, and so on. However, to the best knowledge of the authors, no result 
has been published in the system identification and automatic control field [28]. Thus, this paper focuses on 
presenting the applying possibility of deep learning in system identification areas. 

This study investigates the performance of a one-stage servomechanism system using several 
advanced control techniques. Firstly, identified models using linear least squares, nonlinear least squares, and 
DNN had obtained and selected the best between them. Both the conventional PD and fixed rule base fuzzy 
PD controller had investigated but still, the performance is not acceptable. Therefore, a new variable structure 
(VS) fuzzy PD controller had designed to tune the rule base online using an adaptive mechanism with fast 
calculation to can apply practically. The parameters of the adaptive mechanism had obtained based on the 
selected model reference. The optimal parameters of the proposed control techniques had been determined 
using harmony search (HS) optimization technique. The paper is prepared as follows, firstly, the 
experimental setup is presented. Secondly, system identification is explained. Thirdly, the proposed 
controller techniques are demonstrated. Fourthly, the experimental results are illustrated. Finally, the 
conclusion is discussed. 


2. EXPERIMENTAL SETUP 

This section demonstrates the main parts of the one-stage servomechanism system prototype. Also, 
it shows the open-loop performance of the servomechanism system. Figure 1 illustrates the main components 
of one stage servomechanism experimental setup which consists of a DC motor electro-mechanical module. 
The stroke of the stage ranges from 0 to 9 Inches. 

The DC motor has a nominal speed of 1800 rev/min, and an armature voltage of 90 V dc. The 
optical encoder is an add-on that provides position feedback signal (200 pulses per revolution). Two 
magnetic limit switches detect when the sliding block reaches the start or end position. The DC motor drive 
controls the DC motor electro-mechanical module. This versatile drive also allows an external signal to 
control the motor speed. A data acquisition card (DAQ) NI USB-6009 is connected to the computer that is 
used to perform the control algorithms. The main idea of the program has been designed to make the NI 
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DAQ 6009 generates an analog output signal (-5 to 5 V) to the linear amplifier. Also, the analog output signal 
from the optical encoder has been collected at the same time. The speed of the DC motor will fluctuate when 
the generated signal change continuously. The positive signal will cause the DC motor speed to fluctuate in 
the forward direction, while the DC motor will fluctuate in the reverse direction through the negative voltage 
ranges. The shaft of the optical encoder is coupled with the lead screw shaft to measure the speed and 
position of the stage as in Figure 2. 


One stage table 
- Coupling ftead Scow 
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driver encoder 


Encoder 
pulses 


Optical 
encoder 


and linear 
amplifier 


Push buttons, 
toggle switches 
and lights 


6009 card 


Data collecting 


Figure 1. The one stage servomechanism Figure 2. The block diagram of the experimental setup 
experimental setup servomechanism system 


Figure 3 demonstrates the PRBS output signal from the NI DAQ card. It can be noted that the signal 
has a variable frequency at an acceptable range. Then, this signal will be input to the DC motor driver. Figure 4 
shows the corresponding linear speed of the stage measured by the optical encoder. It is clear that the linear 
speed was fluctuating highly due to the violent excitation of the servomechanism system. Therefore, the stage 
position increases in positive ranges of the input signal entering to DC motor driver while the actual stage 
position decreases in negative ranges of the input signal. The input/output data will be collected and stored in 
an excel sheet file and then this data will be used to develop an identified model for the experimental setup. 
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Figure 3. The random output signal of the NI DAQ Figure 4. The output linear speed of the 
card servomechanism table 


3. SYSTEM IDENTIFICATION 

This section investigates an accurate model for a one-stage servomechanism system based on 
measured input/output data. Three system identification techniques will be used. The first technique uses 
linear least squares. The second technique is the nonlinear least squares. Lastly, the third technique develops 
a DNN-identified model. 


3.1. Linear least squares model 
The general linear transfer function can be expressed as follows: 
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29 — —__ a) 


u(s)  Dy.S+by—1.S"-14---+b9 


In (1), k, bn, ..., bo are the estimated parameters of the approximate transfer function (T.F). The y(s) stands 
for linear stage speed while u(s) is the input voltage to the DC motor drive. The accuracy of the transfer 
function improves significantly when the system degree is increased. However, often there is a restriction 
that increasing order cannot make the model accurate sufficiently [29]. Therefore, it is necessary to explicitly 
add the nonlinearities into the system [30]. 


3.2. Nonlinear least squares model 

As the servomechanism systems suffer from uncertainty and complex nonlinear dynamics, a 
nonlinear ARX (NLARX) model structure has been developed to model such systems [31]. The NLARX 
model consists of several complicated nonlinear functions that can model the backlash and friction in 
servomechanism systems significantly. An NLARX model can be understood as an extension of a linear 
model. A linear SISO ARX model has this structure [32]. 


y(t) + y(t — 1) + azy(t — 2) +--+ anay (t — na) = byu(t) + bou(t — 1) +--+ bapu (t — nb + 1) + e(t) (2) 


Where u is input, y is output, na is the number of past output terms and nb is the number of past input terms. 
This structure can be extended to create a nonlinear form where instead of the weighted sum that represents a 
linear mapping, the NLARX model has a more flexible nonlinear mapping function [33]. 


y(t) = f(y(t — 1), y(t — 2), y(t — 3) .... u(t), u(t — 1), u(t — 2),...) (3) 


Where f is a nonlinear function (to simulate the behavior of friction and backlash that exist in 
servomechanism systems). Inputs to f are model regressors. 


3.3. Deep neural network model 

There are two main types of ANN, the first type is the shallow neural network and the second type is 
the deep neural network. Deep learning has several layers of hidden units and it also permits many more 
factors to be used before over-fitting occurs. Thus, for deep learning, a deep architecture is used. The neural 
network structure in this study contains three layers the first layer is the input layer which receives the input 
PRBS signal to the DC motor driver. The second layer is the hidden layer that contains several hidden 
neurons and receives data from the input layer. The third layer is the output layer which presents the 
corresponding linear stage speed. x; = (x1, X2, X3, ..., Xn)* Input vector applied to the layer, the whole of the 
hidden neuron input ‘j ’ is: 


net}! = ied Wij Xi + or (4) 
With X; = (X4, X2,X3,.,Xn)* (5) 
such as i = 1,2,..., Nh 


x; is the input vector applied to the layer, and w;; is the weights of (i) input neuron connection, and or 
represent the bias of hidden layer neurons. The neurons of the hidden layer can be written as follows: 


yf = fwi Xi tO (6) 


3.4. Candidate identified models validation 

The obtained input/output data was used to develop a set of three candidate-identified models to 
select the best between them. Three models had been obtained, the first identified model is a linear second- 
order system while the second identified model is an NLARX model and lastly, the DNN identified model. 
The linear second-order system has the following specifications. 


y(s) _ 27.05 
u(s)  $2413.91s+78.45 


(7) 


By matching with in (1), the estimated parameters of the linear approximation transfer function are k = 
27.05, n = 2, b, = 1, bn-1 = 13.91 and by = 78.45. On the other hand, the obtained discrete-time ARX 
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model has the following structure. 


A(z)y(t) = B(z)u(t) + e(t) (8) 


Where: 
A(z) = 1 — 0.756427} — 0.09098 z~? 
B(z) = 0.03327 z7! + 0.01742 z7? 

By corresponding with in (2), the polynomial orders are na = 2 and nb = 2. Also, the used nonlinear 
function is a wave net with 25 units. Figure 5 demonstrates the linear speed of the actual experimental setup 
and the candidate models. It is noted that the identified model based on DNN can simulate the behavior of an 
actual experimental setup compared to other identified models. Table 1 demonstrates the mean square error 
of candidate models. It can be noted that the identified model based on DNN has a minimum error compared 
to other identified models. 
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Figure 5. The linear speed of one stage table servomechanism for actual experimental setup and identified 
models 


Table 1. Mean square error of candidate models 


No. System identification method Mean square error 
1 Identified Model (Second Order T.F.) 0.1973 
2 Identified System (Nonlinear ARX model) 0.05912 
3 Identified NARX neural network model 0.00071 


4. VARIABLE STRUCTURE FUZZY PD CONTROL 
4.1. VS fuzzy PD control 

There are two types of adaptive fuzzy controllers. The first type tunes the rules of the fuzzy system 
so, it is called self-organizing controllers or variable structure controllers [34]. In the second type of fuzzy 
controller, the scaling factors are modified online and can be called self-tuning controllers [35]. The previous 
studies showed that the first type of adaptive fuzzy controller is more effective than the second type [36]. 
But, the first type of adaptive fuzzy controller is more complicated than the second type where the first type 
needs building the fuzzy system by the designer without using the fuzzy toolbox for any software program 
[37]. The proposed technique has an adaptive mechanism to tune the centroid of rule base membership online 
based on the optimal model reference adaptive system. Figure 6 illustrates the general structure of variable 
structure fuzzy PD control. 

The variable structure fuzzy PD control input/output relationship can be described as follows: 


u(s) = GU. Au(s) (9) 
Au(s) = 673(x) (10) 
E(s) = GE.e (11) 
CE(s) = GCE. ce(s) (12) 
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GE =1, GCE =, GU =k, (13) 
P 


Adjustment 
Mechanism 


Figure 6. Variable structure fuzzy PD (VSFPD) control 


Where u(s) is the output fuzzy controller, where 87 is the centroid vector of the output membership 
function, 3(x) is the vector of the fuzzy basis function, x € [—1,1], and GU is the scaling factor for the fuzzy 


PD controller. Inputs of the fuzzy controller are scaled using, GCE scaling factors also known as fuzzy gains, 
e(s) is the error between the reference position and the actual position of the stage and ce(s) is the change 
of this error as shown in Figure 7. Here, we introduce a reference model in the structure of the fuzzy 
controller to generate the model error given by: 


€m(S) = y(s) = Ym(S) (14) 


In (7), the y(s) is the actual system output, the y,,(s) stands for the desired performance and e,,,(s) 
considers the difference between the actual process value and the expected value of output. The reference 
model can be a first or second-order system. The model reference transfer function contains the desired 
response of the system such as the desired damping ratio, the desired natural frequency, the desired rise time, 
the desired settling time, and the desired overshoot. If the order of the reference model is a stable first-order 
system as the following transfer function. 


Ym(S) = km 
uc(S) tmst1 


(15) 


Where u, is the reference position, Ym is the output of the reference model, km represents the DC gain of the 
system ratio between the input signal and the steady-state value of output and T, is the time constant that 
measures how quickly a first-order system responds to a unit step input. The MIT rule is the original 
approach to MRAC. The name is derived from the fact that it was developed at the Instrumentation 
Laboratory (now the Draper Laboratory) at MIT. To adjust parameters in such a way that the loss function is 
minimized [38]. 


j(@) = Fe (16) 


To make j small, it is reasonable to change the parameters in the direction of the negative gradient of 
j, that is, 


do aj dem 
a2 _ _y m _ye,, 20m (17) 


where y stands for the adaptation gain while @ is the centroid vector of the output membership function. 
Figure 8 demonstrates the adaptive output membership functions for VSFC. The linguistic labels of the 
outputs are {NB(8,), NM(@2), NS(@3), ZE(@,), PS(@;), PM(8,), PB(@7)}. The output memberships centers 
are not fixed and change continuously through a certain range to prevent the overlapping between the 
memberships. Table 2 summarizes the adaptive rule base of VS fuzzy system. 

From the adaptive output of the membership function can outcome the following equations. 
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0 = [01 82 03 04 0s O6 87] (19) 
The used defuzzification technique is the center of gravity. 


u(nT) = Rj Muay _ 97 x(x) (20) 


rh, u(uy) 


Where u (uj) membership grad of the element uj, u(nT) is the fuzzy control output, n is the number of 
discrete values in the universe of discourse. 


f(E, CE) 


-l poi -0.66 [feat | -0.33 Hi 0 | l 0.33 | I 0.66 | | 1 à 
-1 -0.66 -0.33 0 0.33 066 i nce al ear "a er! S'S Ss 
Figure 7. The input of membership functions Figure 8. The adaptive output membership functions 


(error and change of error) 


Table 2. The adaptive fuzzy rule base 
E/CE NB NM NS ZE PS PM PB 
NB NB(6,) NB(6,) NB(6,) NB(O,) NM(6,) NS(@3) ZE(@,) 
NM NB(6,) NB(@,) NB(@,) NM(6,) NS(@3) ZE(04) PS(@s) 
NS NB(6,) NB(6,) NM(@,) NS(@3) ZE(@,) PS(0;) PM(@,) 
ZE NB(6,) NM(6,) NS(@3) ZE(@,) PS(@s;) PM(@,) PB(67) 
PS NM(6,) NS(@3;) ZE(0,) PS(@;) PM(6,) PB(@,) PB(0;) 
PM NS(6;) ZE(@,) PS(@;) PM(@%) PB(6,;) PB(@,) PB(67) 
PB ZE(0,) _PS(@;) PM(@,) PB(@,;) PB(0,) PB(0,) PB(6;) 


4.2. Harmony search (HS) optimization 

The HS is a phenomenon-mimicking algorithm inspired by the improvisation process of musicians 
[39]. Each controller variable corresponds to each musician [40]. The range of the parameter’s value 
corresponds to the musical instrument’s pitch range [41]. The solution matrix corresponds to the musical 
harmony at a certain time. In this paper, the HS optimization technique was used to obtain the optimal values 
of VS fuzzy PD controller based on MRAS (y is the adaptation gain vector, km is the DC gain, T,,, is the 
desired constant time and the scaling factors of the fuzzy system). The initial population of Harmony 
Memory (HM) is chosen randomly through a preselected range for each controller parameter. HM consists of 
harmony memory solution (HMS) matrix as in (21). 


GUany GEa2) GCEas) Yı aa = 7w Kman Tmaaz 
GU GE GCE y: a K, T, 
(2,1) (2,2) (23) Via, 7 (2, m(2,11) 4m (2,12) 
HM=| . i . eo i o i ; (21) 
GU (yms.1yGE cums,2)GCE ¢ums,3)Y1 amsa) = Y7 ums 10) Km (HMS,11) Tm (HMS,12) 


Similarly, the HM matrix has been established for other proposed controller parameters (classical 
PD control and fixed structure fuzzy control). The objective function corresponds to the audience’s aesthetics 
to be minimized or maximized. The used objective function is shown in (22). 


1 
~ (1-e-P)(Mp+ess)+e~Ê(ts-tr) 


f (22) 
Where M; is the overshoot of system response, e,, is the steady-state error, ts is the settling time and t, is the 


rise time. Also, this objective function can compromise the designer demand by the weighting parameter 
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value (B). The parameter is set larger than 0.7 to reduce overshoot and steady-state error. If this parameter is 
set smaller than 0.7 the rise time and settling time will be reduced. 

The HS optimization program contains several steps. The first step, initialize the HS parameters. In 
the second step, generate random values for the HM matrix and determine the fitness function corresponding 
to each solution vector. The third step, Improvise a new harmony from the HM matrix. The fourth step, 
replace the worst solutions with the best solutions. In the fifth step, repeat step 2 and step 4 until the 
termination criterion is satisfied. Table 3 demonstrates the obtained parameter values using HS optimization 
for the proposed controllers. 


Table 3. The parameter values of the proposed controllers 


Controller Parameter Value 
PD controller ky 26.1677 
ka 1.7924 
Fuzzy PD controller GE 0.087 
GCE 0.0456 
GU 3.01 
VS Fuzzy PD controller km 1.01 
Tm 0.102 
GE 0.092734 
GCE 0.00054574 
GU 40 
Yı 0.04 
Y2 0.045 
Y3 0.056 
Va 0.078 
Ys 0.089 
V6 0.0233 
Y7 0.0123 


5. EXPERIMENTAL RESULTS 

This section demonstrates the motion behavior of one stage servomechanism system using the 
variable structure fuzzy PD control. Also, the performance of the proposed controller technique had been 
investigated by comparing it with the PD control and the fixed structure fuzzy PD control. Two tests had 
been applied to the experimental setup. In the first test, the position command was adjusted at a constant 
value while the second test was executed at the variable position command. 


5.1. Constant position command 

In this test, the position command was adjusted at a constant value (7 inches) and the stage began its 
motion from the position reference (zero inches). Figure 9 shows the position responses of the stage using 
three various control techniques. It can be noted that the variable structure fuzzy PD control has a minimum 
rise time and minimum settling time compared to other control techniques. Also, the variable structure fuzzy 
PD control has no overshoot approximately as demonstrated in Table 4. 

Figure 10 demonstrates the corresponding linear speed responses of the stage through the 
experiment. It can be obvious that the variable structure fuzzy PD control has a high speed at the rise time. 
Also, it is decaying rapidly at moment 2 seconds to compensate for the error in the stage position. Then, the 
stage speed stables at zero approximately. 


8 5 r - ' 
7 =-= HS PD Control 
T4 = = ‘Fuzzy PD Control 
= 6 = ——— VS Fuzzy PD Control 
= <= 
= 3} 
=> = 
Sa = 
Š 3° 
> 2 
$ aE 
a 2 === Position Reference Command S, 
=-= HS PD Control S 
1 — — -Fuzzy PD Control Hw 0 
o ~ VS Fuzzy PD Control 
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Time (sec) : : a Time eee $ 
Figure 9. The stage position responses of proposed Figure 10. The stage speed responses of proposed 
controllers at constant position command controllers at constant position command 
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Figure 11 demonstrates the proposed controller's output through the experiment. The output signal 
in the case of the VS fuzzy PD control referred to the output after the gain GU. It is clear that all controllers 
give the maximum value at the beginning of the test which makes the stage accelerate rapidly to reach the 
desired position. Figure 12 illustrates the online change in the output membership centers through the 
constant position command test. It is noted that the adaptive mechanism tunes (81, 02, 03, and 04) where the 
fuzzy output value locates through these memberships so, the other membership centers (85, 06, and 07) are 
not changed. 


Table 4. The performance at constant position command test 


Controller type Rise time (sec) Settling time (sec) Overshoot (%) 
The HS PD control 2.7123 3.5678 0.2123 
The fuzzy PD control 1.6512 2.6565 2.0702 
The VS fuzzy PD control 1.2344 1.8343 0.1778 
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38 ` Ee 2 ee 
` 1 H a15 
2 ` e—a ee 
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Figure 11. The stage position responses of proposed Figure 12. The change in membership centers during 
controllers at constant position command the constant position command test 


5.2. Variable position command 

The position command in this test changes continuously to investigate the robustness of the 
proposed controller technique and its ability to track a complex position command profile. Figure 13 
illustrates the position responses of the stage by several control techniques and under variable position 
command (half-sine wave). It is obvious that the variable structure fuzzy PD control can follow accurately 
the complicated trajectory compared to other control techniques. Also, the fixed structure fuzzy PD control 
has acceptable tracking accuracy but it has a high error at the beginning of tracking while the harmony search 
(HS) PD control has a high deviation from the position command profile. Moreover, the variable structure 
fuzzy PD control has the smallest mean square error throughout the experiment as demonstrated in Table 5. 
Figure 14 demonstrates the corresponding linear speed responses of the stage for control techniques. It can be 
noted that the stage speed through the variable structure fuzzy PD control has a high value in the first seconds 
of the experiment and then the stage speed decreases gradually to track the position command profile. 
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=~ = — ‘Fuzzy PD Control 
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Figure. 13. The stage position responses of each Figure 14. The stage speed responses of control 
control technique at variable techniques at variable position command 
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Table 5. The mean square error of proposed controllers at variable position reference command 


Controller Type Mean square error (inch) 
The PD control 0.912 
The fuzzy PD control 0.7065 
The VS fuzzy PD control 0.412 


Figure 15 shows the controller's output through the variable position command test. At the start of 
the experiment, the controllers have a high output value. These signals decrease gradually with different 
behaviors. The controller's signal declines until the polarity of signals change which leads to the reverse of 
the stage motion direction. Figure 16 shows the online change in the output membership centers during the 
variable position command test. It is obvious that the adaptive mechanism changes (01, 02, 03, and 04) 
where the fuzzy output value belongs to these memberships so, the other membership centers (05, 06, and 
87) do not change until the end of the test. 


The change in membership centroid 
o 


Controller Output (volt) 


o os 141 145 2 25 3 35 4 45 5 Time (sec) 
Time (sec) 


Figure 16. The change in membership centers during Figure 15. The control techniques output at variable 
the variable position command test position command 


6. CONCLUSION 

A new variable structure technique has been developed to tune the rule base of fuzzy logic control 
online based on optimal model reference adaptive system (MRAS). This work investigates the robustness of 
the proposed technique by applying it to one stage servomechanism system. The controller aims to track 
accurately a preselected position reference trajectory in presence of the friction and backlash problems. Also, 
the performance of the proposed controller technique has been compared to the classical PD control and the 
fixed structure fuzzy PD control to ensure robustness. There are two tests that have been implemented to 
investigate each control technique. The first test adjusts the position reference at a constant value while the 
second test adjusts the position reference to change continuously with time. The experimental results 
demonstrate that the variable structure fuzzy PD control based on an optimal model reference adaptive 
system can eliminate the tracking error quickly compared to other control techniques which will save the 
machining time cycles. 
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